Pure-radiation gravitational fields with a simple twist and a Killing vector 
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Pure-radiation solutions are found exploiting the analogy with the Euler-Darboux equation for 
aligned colliding plane waves and the Euler-Tricomi eqution in hydrodynamics of two-dimensional 
flow. They do not depend on one of the spacelike coordinates and comprise the Hauser solution as a 
special subcase. 
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2 ■ I. INTRODUCTION 

fS| [ There exist many papers dealing with algebraically special, expanding and twisting pure-radiation solutions of the 
5^ Einstein equations. An extensive bibliography up to 1980 exists in Ref. [1]. Further results on pure-radiation fields 
can be found in Refs. [2-7]. The standard form of the metric is [1] 

f-^ : ds'^ = - 2n [dr + WdC + WdC + HQ] , (1) 



(N 



0^ 



du + LdC + LdC. 



Here r is the coordinate along the null congruence of geodesies, u is the retarded time, while C span a two- 
' dimensional surface. The metric components are determined by the r-independent real functions P, m, M and the 
1 complex function L 



2iE = p2 (^Q]^ _ ^2.) , (2) 



O 

o 

I ' r + 

5-1 . 

W = p-'^Lu + id^, (4) 

• 1— I 

X, 

H = -r {In P)^-{mr + MJ:)pp+-, (5) 

K = 2P^Re[d{dlnP-L,,)], (6) 
where 8 = 8^ — Ldu and E is the twist. The functions mentioned above satisfy the system of equations 

[d - SLu) (m + iM) = 0, (7) 

P-^M = ImddBBV, (8) 

= -2P3 [p-3 + iM)]^ + 2P3 [ddddV)^ - 2P^ (ddV)^ {9dV)^ , (9) 



where Vu = P , n is the energy density of pure radiation and the Newton constant is set to 1. Equations (7)- (9) are 
in fact Eqs. (26.32) and (26.33) from Ref. [1]. 

It has been noticed in different contexts that the condition M = simplifies the equations [4-6,8,9]. In a previous 
paper [10] we explored this condition with the help of Stephani's method [2] for some algebraically special, axisym- 
metric, expanding and twisting gravitational fields. They depend on u and cr = CC intrinsically i.e. the w-dependence 
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cannot be taken away by applying the gauge transformations (25.27) from Rcf. [1]. In the present paper we find 
solutions for fields which depend on u and the real part of^, x = "^(C + C) ^iid have the simplest possible twist. 
Many of the well known algebraically special exact solutions possess this kind of symmetry [1,3,4]. Among them 
arc the only known seed vacuum Robinson- Trautman solution of type III. the static C-mctric and the only known 
expanding vacuum solution of type N with twist, namely the Hauser solution, given respectively by Eqs. (24.15), 
(24.23) and (25.71) from Ref. [1]. 

In Sec. II Eqs. (7)-(9) arc reformulated in terms of an invariant potential which leads to the Lu = gauge. In Sec. 
Ill solutions with separated variables ?i and x are found. In Sec. IV the main Eq. (8) for simplest twist is shown 
to be equivalent to a case of the Eulcr-Darboux equation, somewhat different from the equation derived in Ref. [10]. 
On its turn it is a complex version of the central equation in the theory of aligned colliding plane waves (CPW). 
Three new solutions are found. In Sec. V a homogenous hypergeometric solution is derived, exploiting the similarity 
between Eq. (8) and the Eulcr-Tricomi equation. It contains as a special case the Hauser solution of type N. Section 
VI contains some conclusions. 

II. FIELD EQUATIONS IN THE Lu = GAUGE 

Following [2,10] we introduce the invariant complex potential ^ which solves Eq. (7): 

m + iM = (f)l, (10) 

\/2(j)u ^ ^ 

When M = we can apply the transformation 

u' = f{u,x), (12) 

(to + iM)' = (to + iM) (13) 
to make m a positive or a negative constant TOq, so that 

cP = m'J^[u + iq{x)], (14) 

L = -^fe (15) 

with real q. Obviously i„ = 0. This gauge differs from the usual gauge P„ = but is very suitable when the 
Newman-Unti-Tamburino (NUT) parameter M vanishes. Equations (8) and (9) simplify 

ddddV = BdddV, (16) 

= 6tooP"^P„ + 2P^ddddP - 2P'^ddPddP (17) 

with 9 = -1= (dx — iQxdu)- The second equation is in fact an inequality. When Pu ^ ri^ can be made positive by 
the choice of too at least for some region of spacetime [1,4,8]. The expressions for the metric components simpliiy 
too, e.g., the gauge invariants S and K become 

S = ^fcP', (18) 

K P"^ {dd + dd)\nP. (19) 

When TO + iM = (Pctrov types III and A^) Eq. (7) is an identity but still a potential (/> may be introduced with 
the property dcp = and the subclass of solutions satisfying Eq. (14) (with toq = 1 ) can be studied. One should put 
mo = in all other equations. 
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The main equation (16), which is of fourth order with respect to V, becomes in both cases a hnear second order 
equation for P. Equation (18) shows that q must be at least quadratic in x for a non-trivial twist. Let us choose the 
simplest possibility: q = L = -^x. Then equations (16) and (17) read 

x'Puu + Pxx = 0, (20) 



2n2 = 12moP-'P„ - 3p3p„„ - ^x^P^P^ - P^ (P„ + 2a;P„^)^ . (21) 

The last two terms in Eq. (21) are definitely negative, so the first term must be necessarily positive for a type II 
solution and the second term must be positive for a type III solution. 



III. SOLUTION WITH SEPARATED VARIABLES 

Suppose that P = F {x)G {u). Equation (20) splits into two parts 

Gun = cG, (22) 

F,, + cx^F = 0, (23) 
where c is an arbitrary constant. There are three types of solutions, depending on the sign of c. If c = we have 

P = ux, (24) 

= ^ - '-u'x\ (25) 

Type III solutions have negative energy while type II solutions have negative energy density for x ^ oo. If c > 0, P 
contains Bessel functions, e.g., 



(26) 



= -2A^/Zmo - 2c (3 + 4cx^) F^e'^^^" - 2cP'^ {F + 2xF^f e'^^". (27) 

The first term can be made positive when mo < 0, but the second has a negative pole when x ^ oo because 
{xF) ~ x"^. There are no solutions with everywhere positive n^. The third case c < gives a generic solution with 
modified Bessel functions like 



P = ^JxKi/i -^-^r— sin yf^cu, (28) 



= 24V^mo cot {^f^cu) - 2c (3 + 4ca;'') F'^ sin'' ^f^cu + 2cP^ (P + 2xF^f cos^ (29) 

The first term changes sign and has poles in u, thus type 11 solutions are unphysical. Type 111 solutions have negative 
'n? when e.g. u — tt/-/^, although the second term doesn't have poles in u. 

Equations (24) and (26) have been obtained by another method in a more general form in Ref. [4] but the energy 
density has not been discussed in details. 

If we separate the variables like P = F (x) + G (u) then 

cu 

G= — — + C1U + C2, (30) 



P = ^ + C3 + C4, (31) 

where c^ are constants. Numerous negative terms arise and the region of positivity of is rather complicated due to 
the many arbitrary constants. A more detailed discussion in a general setting of this type of separation of variables 
may be found in Ref. [11]. 
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IV. REDUCTION TO THE EULER-DARBOUX EQUATION 



Equation (20) with replaced by where / is an integer has been studied in the past [12,13]. Green functions for 
different boundary problems have been found. The basis of these results is the change of variables which transforms 
Eq. (20) into the Euler-Darboux equation 



1 



<^ + w 



(P^ + P^) = 0, 



(32) 



(33) 



where ui = — (^.This equation is similar to Eq. (24) from Ref. [10] and the main equation for aligned colliding plane 
waves [14] - one needs only to replace the multiplier 4 by 2. This, however, is not a trivial change and can't be 
achieved by scaling the variables or P. Nevertheless, one may use the numerous techniques developed in the search 
for CPW solutions and applied in Ref. [10] for axisymmetric fields. We shall find analogues of these solutions. As 
pointed out in Ref. [10] the reality of P must be ensured because the variables y, co are complex. 
The simplest solution 



is a time-independent vacuum solution of Kerr-Scliild type. 

A solution with separated variables P = F {tp) G (to) exists. Replacement in Eq. (32) yields 

P = A[{a + ^) {a + ^)]-'/\ 
where A and a are constants. A is ignorable and a can be hidden in u to obtain a real solution: 

P = 



(34) 



(35) 



(36) 



B = «2 + _a;4. 

4 



(37) 



This solution is analogous to the solution given by Eqs. (27) and (28) from Ref. [10] and has a number of nice features. 
The energy density is 



n 



53/2 



128B5 



4x4 {6u^-x^Y + u^ {9x^-4u^y 



(38) 



In the following we suppose that the retarded time satisfies the condition u > uq > for some constant uq. Equation 
(38) is regular in x unlike many other solutions, plagued by singular pipes for x = or a; = ±00 [1]. When mo < 
the first term dominates over the others if | mo | is big enough and consequently is positive. Unfortunately, type 
III solutions arc not with positive for any x because the second term changes sign. The gauge invariants (18) and 
(19) are regular in x and vanish when u 00: 



(39) 



K = 



453/2 • 

The same is true for the Weyl scalars [1,15] with leading terms given by 

*2 = mop^ 

*4 = pP^Iu + O , 



(40) 

(41) 
(42) 
(43) 
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/ = p-^ddP. (44) 

An exception is ^'2 which approaches —nio/r^ as u ^ 00. 

Let us present next an analogue of the cosh~^ solution found in Ref. [10]. We substitute the ansatz P = P (a) , 

a=^^^ = ^, (45) 

into Eq. (32). The result is an elliptic integral of the first kind F {ijj, k) : 

P = V2F (^V, ^) , (46) 

1I) = arccos ( 1 + ) ~ ^^"^ , (47) 

Pa = {l + a'y'"\ (48) 
The function P is bounded: < P < V2F (f , 75)- 

The solution possesses regular characteristics like the previous 

one: 

2 _ 12mo I X I 18m I a; I P3 9a;2p2 
" ~ P(4i3)^/^ ^ (4P)^/^ ~ ^S^' ^ ^ 

E = \P\ (50) 

^ = -;^- (51) 

All terms in Eq. (49) are regular. The first term is positive for mo > and type II solutions with positive exist for 
big enough niQ. The second term is positive too but does not always dominate over the third one. The Weyl scalars 
are also regular and ^4 vanish when u ^ 00. 

Let us transform now the Euler-Darboux equation (32) into its canonical form. Introducing the new variables 
T = x^, A = 2u we obtain 

Prr + ^Pr + P\X = (^- (52) 

It 

This is an analogue of Eq. (25) from Ref. [10] and its solution is given by the Bessel functions from Sec. III. We can 
go further, utilizing the coordinates for the first Yurtsever solution [10,14,16]: 

T = vsm.r}, (53) 
A = vcosr], (54) 
= 4B, (55) 

cos77=(l + ^) . (56) 



Then Eq. (52) becomes 



P.. + ^Pr>r> + \ i^P. + ;^ cot 77 P^ ) = 0. (57) 
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It has a separable solution of the kind P = v^'Y (77). Y, instead of being a Legendre function of the first or second 
kind, satisfies the equation 



(1 - w^) Y^^ - ^wY^ + l(^l+l^Y = 0, 



where w = cos rj. Its solution is a hypergeometric function and 



where e + cr = — |, sa = I (^l + ^) and 



P={4By/^F[s,a,~,X], 



X 



P is real because < X < 1. The hypergeometric function is reducible to a Legendre function: 

„4 \ 7/8 



(58) 



(59) 



(60) 



(61) 



If Z = — 5, either e or a vanishes and formula (59) degenerates to the rational function given by Eq. (36). 



V. HYDRODYNAMICAL ANALOGY 



Equation (20) has certain similarities with the Euler-Tricomi equation 



uu ^xx 



^x 



0. 



(62) 



It appears in hydrodynamics in the study of a two-dimensional flow of compressible fluid with velocity near the 

velocity of sound [17]. It is a limiting case of the more complex Chapligin equation which has integrals among 
the hypergeometric functions [17,18]. The Euler-Tricomi equation is invariant under the transformations —> cu^, 
cx^ which leads to a homogeneous hypergeometric solution. It is discussed at length in Ref.[17]. In our case 
Eq.(20) is invariant under u (?u, x ^ cx and we can try a homogenous solution 

P = u'^F {z) , (63) 



z = — 



where k is the degree of homogeneity. Then Eq. (20) becomes 

ziz-l)F,,+ f^-k\z-^ 



Once again this is a hypergeometric equation and one of its fundamental solutions leads to 



P = u''F[ -- 



k l-k 3 



2' 2 '4' 

In fact, the hypergeometric function in Eq.(66) degenerates to a Legendre function for any k 

[(i-.)-v^l 



P = 2-V4r ( ^ ) (-^)i/8 (1 _ pV4 



It becomes a rational function in some cases. Thus if fc = ^ 

1/4 



P=(|)-[i+(i-.)V^; 



fe-3/4 



1/4 



(64) 



(65) 



(66) 



(67) 



(68) 
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and if fc = — I 



1/4 



Exploiting only Eq. (64) one finds the following identities 



.1/2 



1/4 



Puu — Pl^O ' 



(69) 

(70) 
(71) 
(72) 



Pi = 4k{k-l)P- (4fc - 3) xPa: 



(73) 



They hold for any solution of Eq. (65). Plugging Eqs. (70)-(73) into Eq. (21) one can study the properties of the 
energy density. It simplifies drastically when A; = |: 



(74) 



= moP-^Pu 



mo 
4u 



V8' 8' 4' 
\ 8' 8' 4' ^/ 



(75) 



When <| z |< 1 it can be checked with MAPLEV that the r.h.s. of Eq. (75) is positive for mo > 0. In the region 
1 < I -2 I < oo an analytic continuation of the hypergeometric functions should be performed. Thus 



3 1 3 



13 3 1 



(76) 



r(!)ra) 
ra)r(l)' 



A2 = 



r(l)r(-^) 
r(-f)r(|)' 



and similarly for the other hypergeometric function in Eq. (75) . A check with MAPLEV confirms again the positivity 
of n^. The energy density is regular in u and z. When mo = a vacuum solution is obtained. It has ^'2 = and then 
the higher terms in Eq. (42) vanish [1] . Prom Eq. (44) it follows that 

(77) 



4 {x^ - 2iu) 

and dl = 0, J„ 7^ 0. Therefore = 0, while ^'4 ^ and the solution is of type A''. In fact, this is the Hauser solution 

[1,19] in the Lu ~ gauge. 

Let us show this in detail. We have used the transformation (13) to bring (j) to the simple form (14) (with mo = 1 
since the field is of type N) and then have dropped the primes. Restoring them, Eqs. (12)-(14) show that 

f = u' = Cx^u, (78) 

where C is an arbitrary constant which we fix to C = 1/2. Under transformation (12) P and L change as 
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(79) 



L' = fuL - -^U, (80) 

(see Eq. (25.27) from Ref. [1]). P' is given by Eq. (74) while L' = ix/V^. Then the original P and L are given by 

P = x'^/^F (u) , (81) 

L = ^{u + i), (82) 

X 

F in) = Uf (-1-11, -uA + \a,uF ill l-uA . (83) 



2 V 8 8' 2' y 4 V8 8' 2 

Wc have used formula (76) to derive the expression for F (u). Equations (81)-(83) give exactly the Hauser solution 
[19] as written out in Refs. [1] and [20]. Equation (83) is a specific linear combination of the even and odd solutions 
given in Ref. [20] which fixes the one-parameter freedom of the Hauser solution. This ends the proof of our assertion. 

VI. CONCLUSION 

We have shown that when the NUT parameter M vanishes and the gauge L„ = is used, the main equation (20) for 
expanding pure radiation fields with a simple twist and a special symmetry becomes a tractable second order linear 
equation for P. It is reducible to a case of the Euler-Darboux equation, somewhat different from the central equation 
in the theory of aligned colliding plane waves. We have found regular solutions, separating the variables in different 
coordinate systems. In some cases the regions of positivity of the energy density were investigated. Another analogy 
with the Euler-Tricomi equation, appearing in the hydrodynamics of two-dimensional fluid flow, has been exploited 
to find homogenous solutions. Interestingly enough, the Hauser vacuum solution of type N is an exceptional member 
of the family of type II solutions with degree of homogeneity | . 
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